Which Chiral Symmetry is Restored in High Temperature QCD? 
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Sigma models for the high temperature phase transition in quantum chromodynamics (QCD) 
suggest that at high temperature the SU(Nf) x SU(Nf) chiral symmetry becomes exact, but the 
anomalous axial (7(1) symmetry need not be restored. In numerical lattice simulations, traditional 
methods for detecting symmetry restoration have sought multiplets in the screening mass spectrum. 
However, these methods were imprecise and the results, so far, incomplete. With improved statistics 
| and methodology, we are now able to offer evidence for a restoration of the 5(7(2) x 5(7(2) chiral 

^ . symmetry just above the crossover, but not of the axial (7(1) chiral symmetry. 



^T) , I. INTRODUCTION 

o , 

A high temperature phase transition from a deconfined quark plasma to a confined phase is thought to have occurred 
as the early Universe cooled. This phenomenon is under investigation in high energy heavy-ion collisions. Through 
^\ • numerical simulations of quantum chromodynamics (QCD) we hope to gain an understanding of the qualitative and 
quantitative characteristics of this phase transition. The phase transition (perhaps only a crossover at physical quark 
i masses) is associated with the spontaneous breaking of the chiral symmetry and formation of chiral condensates. 
Sigma models suggest that in the limit of zero up and down quark masses, the SU(2) x 5(7(2) chiral symmetry is 
exact in the high temperature phase Q , and a phase transition separates it from a cold phase in which this symmetry 
is spontaneously broken. The gauge anomaly, present at low temperature, may persist at high temperature, however, 
r* ■ breaking the (7(1) axial symmetry at all temperatures. 

Early efforts to detect symmetry restoration looked for chiral multiplets in the screening mass spectrum For 
example, the following channels are related according to the indicated symmetries: 

<- 517(2) x 5(7(2) -> 

T 

U(1)a fa 

I V ao 



The screening mass spectrum is found from the space-like hadron propagators. The restoration of the 5(7 (2) x 5(7(2) 
symmetry requires a degeneracy between the lowest pion screening mass and that of its chiral partner, the J p = + , 
I = /o meson (also known as the a). The determination of the fo screening mass through numerical simulation is 
complicated by the presence of quark- line disconnected graphs. Computing them requires an expensive determination 
of the quark propagator from multiple origins. In early simulations, therefore, it was common to keep only connected 
graphs. This practice, applied to the fo, results instead in a determination of the screening mass for the J p = + 
I = 1 ao meson (also known as the S) S. This meson is the axial (7(1) chiral partner of the pion. Thus a degeneracy 
in the it and ao screening masses would imply a suppression of the gauge anomaly and a partial restoration of the 
axial (7(1) symmetry, but does not test restoration of the 5(7(2) x 5(7(2) symmetry. 

New simulations with large data samples make it possible to revisit the question of which symmetry is restored 
Further statistical improvement can be obtained by studying the susceptibilities related to the propagators, 
rather than just the screening masses: for example, from the pion susceptibility 
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X, = / d 4 r (7r(0)7r(r)> 



(1) 



and the related susceptibilities, Xfo an d Xa , we can define two order parameters 

Xsu(2)xsu(2) =Xn~ Xfo and XU(1) =Xtt- Xa - (2) 

Restoration of either symmetry requires that the corresponding order parameter vanish. 

We use the staggered fermion scheme. This scheme breaks all but one generator of chiral SU(A) x SU(4). The full 
symmetry is expected to be recovered in the continuum limit. The one surviving generator, however, can be used to 
explore symmetry restoration at the phase transition at nonzero lattice spacing. The staggered fermion treatment 
of the axial U(l) symmetry is less satisfactory. That symmetry, formulated in the conventional manner, is broken 
explicitly on the lattice. It, too, is expected to be recovered in the continuum limit. Since our analysis treats only 
one lattice spacing, namely a ps 1/(6T c ), further study will be required to distinguish between effects of the lattice 
approximation and continuum effects of the gauge anomalies. 

A preliminary report of our results was presented at Lattice '96 || . A number of other groups have also taken up 
this question and have also reported preliminary results |(3|-[|] . 



II. FORMALISM AND COMPUTATION 



We simulate the iV^-flavor staggered fermion action with the standard partition function at temperature T on a 
hypercubic Euclidean lattice with spacing a, quark matrix M(U, m q ), quark mass m q , and gauge link matrices U 0: 

z = e ~VF(T,arn q )/T = f [ dU ] e ^ S g (U)][det M (U, m q )] N f ^ . (3) 



As is well known, the fermion determinant can be expressed as det M(U,m q ) = det[D 2 + (2am q ) 2 ], where the latter 
determinant is taken on the even lattice sites only and D 2 is the square of the fermion hopping matrix. Thus the free 
energy is manifestly even in the quark mass. 

We will be concerned with a variety of susceptibilities related to the singlet chiral order parameter, 

</o> = (W) = dF(T, m q )/dm q = TN f a/2V <Tr M~ l ) , (4) 

where the expectation values are defined on the ensemble (||) . The associated susceptibility is 

Xfo = 9 (f ) /0m q = J d 4 x[{fo(0)fo( X )) - </ (0)) 2 ] = Xconn + Xdisc (5) 

The quark-line connected and disconnected contributions are 

Xco„n =TA^ /a 2 /^(TrA/- 2 ) and X disc = T/V [<(oJV>/2Tt A/" 1 ) 2 ) - (aN f /2Ti Af" 1 )'] (6) 

It can be seen from this result that the disconnected contribution to the susceptibility is just proportional to the 
"configuration variance" of (/ ), that is Xdisc = V/T (/ 2 ) — (/o) 2 • 

All of our simulations are carried out with two dynamical (sea) quark flavors. However, in measuring susceptibilities, 
we can adjust the valence flavor number to suit the observable. If we stick with only the four flavors forced upon 
us by fermion doubling in the staggered fermion scheme, all isospin components of the ao meson are generated by a 
nonlocal fermion bilinear fic|| . However, at the expense of increasing the flavor degeneracy to eight, we can create an 
ao analog from a diagonal fermion bilinear operator. In any case all such ao components are expected to be degenerate 
in the continuum limit and any of them can be used to test symmetry restoration. The susceptibility of the diagonal 
ao operator is exactly the connected part of the /o susceptibility: 

Xa — Xconn • (7) 
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FIG. 1. Phase diagram for the standard SU (3) 
Wilson gauge plus two-flavor staggered fermion 
action showing the approximate Nt = 6 crossover 
location (crosses and burst) as a function of gauge 
coupling 6/g 2 and quark mass am q . Data sample 
points are indicated by octagons. 
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FIG. 2. Chiral order parameters extrapolated in 
quark mass squared. 



We measure this susceptibility directly from the connected part of the fo correlator: Xconn = / d 4 x (/o(0)/o(r)) | conn , 
while Chandrasekharan and Christ measure it by taking the derivative of (fo) with respect to the valence quark mass 
H . Finally, a well-known Ward identity relates the pion susceptibility to the chiral order parameter O] : 



N f Ta 2 /V(Tr{M^M)- 1 ) = (f ) /(2m ? ) 
In practice we measure the order parameters (|^) through 

(fo) /(2m„) 



XS , (7(2)x5i7(2) 



Xdis 



and 



Xc/(i) = </o) /(2m g ) - Xc 



(8) 



(9) 



The simulation consisted of a subset of configurations generated in an extensive study of the equation of state for 
N t = 6 and N f = 2 at 6/# 2 = 5.45 and quark masses am q = 0.0075, 0.01, 0.0125, 0.015, 0.02, and 0.025 |||. This 
parameter range lies in the high temperature phase slightly above the phase transition, as illustrated in Fig. [I], and 
was selected to permit an extrapolation of the measured quantities to zero quark mass in the high temperature phase. 
The simulation sample at each mass covered a molecular dynamics time span of at least 2000 time units with the first 
400 omitted. Measurements were taken at intervals of at most 50 time units. The chiral order parameter (fo) = (^ip) 
was measured using the random source method [ fl2[ with 33 random sources. These measurements, with care taken to 
avoid biases inherent in the noisy source technique, in turn, provided an estimate of Xdisc through the configuration 
variance. 



III. RESULTS AND CONCLUSIONS 

Results are shown in Fig. and table |. We have indicated a linear extrapolation in (am q ) 2 . Because they are 
closer to the crossover (Fig. |lj), where curvature may be expected, we chose to exclude the two highest mass points 
from the fit. The zero mass intercepts are 

Xsc/(2)xsc/(2) =0.04(31) and X u(i) = 0.75(22) (10) 

with x 2 /df = 2.6/2 and 2.5/2 respectively. Fits to all points gave Xsu(2)xsu(2) — —0.33(20) with x 2 /4f = 5.6/4 and 
Xu(i) = 0.81(11) with 2.7/4. 

It is surprising that a fit of the same points to an expression linear in am q gives a result consistent with a zero 
intercept for both order parameters: Xsu(2)xsu(2) = 0.15(38) with x 2 /rf/ = 1.8/2 and Xu(i) = 0.40(56) with 2.4/2. So 
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which fit is correct? As we have emphasized, the free energy is rigorously even in the quark mass. In consequence the 
order parameters are also even. Thus if the free energy is analytic at zero quark mass, a quadratic fit is required. Now 
some gauge field configurations give rise to fermion zero modes or near-zero modes. In a two-flavor simulation, those 
modes contribute terms in [(am g ) 2 ]' v "' 4 = \am q \ to the free energy - terms linear but nonanalytic. Such behavior, 
if not suppressed by a vanishing probability for encountering zero modes, would imply a phase transition or infrared 
singularity at zero quark mass. However, measurements of screening masses for T > T c give no indication of infrared 
singularities for small am q . A phase transition at zero quark mass for T > T c is likewise unexpected in sigma models. 

In conclusion, our results are consistent with the sigma model scenario: a restoration of SU(2) x SU(2) but not 
of U(1)a (approximately 3cr). Whether the apparent breaking of the axial U(l) symmetry is a lattice artifact or a 
consequence of the anomaly remains to be established by future measurements at smaller lattice spacing and with 
improved actions. 
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TABLE I. Susceptibilities and order parameters in lattice units. 



am, 




Xconn 


Xu(i) 


Xdisc 


XSU(2)xSU(2) 


0.0075 


0.0446(12) 


5.21(17) 


0.74(23) 


0.89(21) 


0.15(31) 


0.01 


0.0599(16) 


4.61(9) 


1.38(18) 


0.91(12) 


-0.47(22) 


0.0125 


0.0724(16) 


4.35(9) 


1.44(16) 


1.25(18) 


-0.19(22) 


0.015 


0.0885(15) 


4.21(7) 


1.69(12) 


1.12(20) 


-0.57(23) 


0.02 


0.121(5) 


3.59(14) 


2.5(3) 


3.1(1.0) 


0.7(1.1) 


0.025 


0.157(3) 


3.04(8) 


3.23(14) 


3.3(5) 


0.1(6) 
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